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Explicit results at tree and one-loop order, obtained by the HELAC/PHEGAS package that is based on the Dyson- 
Schwinger recursive equations approach, are briefly presented. 



1. Introduction 

Recursive equations to compute scattering ma- 
trix elements have been used extensively over the 
last years in order to obtain results for multi- 
leg amplitudes. Their history started essentially 
with the work of Berends and Giele , who were 
able to prove the conjectured simple all-n form 
of Parke and Taylor [5] for the MHV amplitudes 
in QCD. The recognition of their usefulness has 
been expanded recently by the discovery of a new 
class of recursive equations, by Britto, Cachazo 
and Feng g] and Witten g]. 

In this paper we are considering the Dyson- 
Schwingcr (DS) recursive approach |5lril7l8j . and 
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show how this can be used as a general frame- 
work for scattering elements computation. We 
also present selected results for processes at tree 
order and at the one-loop level, obtained with the 
HELAC/PHEGAS Jj package, which is an implemen- 
tation of the DS method. 

2. The Dyson-Schwinger approach 

The traditional representation of the scattering 
amplitude in terms of Feynman graphs results to 
a computational cost that grows like the number 
of those graphs, therefore as n\ (at tree order), 
where n is the number of particles involved in the 
scattering process. 

An alternative 1 to the Feynman graph repre- 
sentation is provided by the Dyson-Schwinger ap- 
proach [Jj. Dyson-Schwinger equations express 
recursively the n-point Green's functions in terms 
of the 1—, 2—, . . . , (n — l)-point functions. In the 
framework of a theory with three- and four-point 
vertices the DS equations are rather simple and 

1 For other alternatives see 1101111 . 
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their diagrammatic representation is given below, 
for 1 -> n |12I13I14I15I16I17| amplitude: 




Omitting the contribution of the second line in 
the above formula is equivalent to restrict our- 
selves at tree order. In order to get an idea of the 
actual mathematical form of these equations, let 
as consider the simplest case where we are inter- 
ested to " count graphs" , so by dropping all propa- 
gators, couplings, wave-functions, etc, we end up 
with the following equation: 

a(n) = 5 n ,i + -z: — 7—ra(ni)a{n 2 )S ni+n2 n 
21 * — ' n\\n2\ 

+ ¥i ^2 — T~i — :a,(n 1 )a(n 2 )a(n 3 )S ni+n2+n3 , n 
61 z — ' niin2 ] .n3 ] . 

with the initial condition a(0) = 0; a(n) is noth- 
ing more than the number of Feynman graphs, 
contributing to the 1 — > n matrix element. 

The computational cost of DS equations grows 
like ~ 3™, which essentially counts the steps used 
to solve the recursive equations. Obviously for 
large n there is a tremendous saving of compu- 
tational time, compared to the n! growth of the 
Feynman graph approach. 

2.1. Color representation 

Color representation or color decomposition of 
the amplitude is a major issue when dealing with 
multi-parton processes. Let us consider n-gluon 
scattering with external momenta {pi}\, helici- 
ties {£i}i anQl colors {ai}™ of gluons i = l,...,n. 
As is well known the total amplitude can be ex- 
pressed as a sum of single trace terms JB] : 

M({p i y?,{e i y?,{a i }?) = 2ig n - 2 
Tr(t a H a °i™ . . .t°-*<»)).Aj({p,}?,{e,}?) 

/eP(2,..,n) 



where 07(2 : n) represent the 7-th permutation 
of the set {2, . . . , n} and Tr{t ai t a "i^ . . . t a "i^ ) 
represents a trace of generators of the SU(N C ) 
gauge group in the fundamental representation. 
For processes involving quarks a similar but much 
more cumbersome expression can be derived |18| . 

One of the most interesting aspects of this de- 
composition is the fact that the Ai({pi}i, {&t}" ) 
functions (called dual, partial or color-ordered 
amplitudes), which contain all the kinematic in- 
formation, depend on the permutation and are 
gauge invariant and cyclically symmetric in the 
momenta and hclicitics of gluons. The color or- 
dered amplitudes are simpler than the full ampli- 
tude because they only receive contributions from 
diagrams with a particular cyclic ordering of the 
external gluons (planar graphs). 

Of course to get the full amplitude one has to 
square the matrix element, 

J2 IMte}M^}?,R}i)l 2 

e ij 

where the (n — 1)! x (n — 1)! dimensional color 
matrix can be written in the most general form 
as follows: 

Cu= Tr{t a H a °iV> . . .t a "i^)Tr{I <-► J)* 

l..JV c 

(1) 

There exists a much simpler approach, in fact far 
superior from the point of view of an automatized 
numerical calculation, where the matrix element 
is represented as follows |19I7I2(J|5] . 

<»<}?) = 

2ig n ~ 2 ]T Dj A I {{ Pl } n ll {e l } n 1 ) 

I=P(2,...,n) 

with Ci , ai the color and anticolor indices for 
each external particle, i.e. (c, 0) for quarks, (0, a) 
for antiquarks, (c, a) for gluons and (0, 0) for non- 
colored particles, and 
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or in a more abstract notation 

Dl = §l,a-j{l)^2,ai(2) ■ ■ ■ bn,cn(n) 

where 07(1 : n) represent the 7-th permutation 
of the set {1,2, ... ,n}. The sequence of num- 
bers i, <Ji(i), i = 1 . . . n, is identified as a color- 
connection configuration, describing the way the 
color connection is structured. In that sense, no 
explicit reference to 'real' color indices is made. 
Finally the color matrix takes a very simple form, 

Cij = N™ {a '-' 7j) (2) 

where 1 < m(ai,aj) < n counts how many com- 
mon cycles the permutations 07 and trj have. For 
a detailed description, see 0. 

Recursive equations can be written both for the 
full amplitude, M., and for the color ordered, A. 
In the latter case the DS equations are identical 
to the Berends-Giele ones. 

For numerical applications the computation of 
the color ordered amplitudes suffers from the n! 
growth related to the number of color-flow or 
color-connection configurations. In such cases 
it is preferable to write down DS equations for 
the full amplitude M({pi}i, {ei}™, {ci, ai}") and 
then perform the incoherent sum 

Ci,ai=1...3 

by Monte-Carlo methods. We have recently 
extended HELAC so that a Monte-Carlo over 
'real' colors, or color-configurations can be per- 
formed A color configuration is identified by 
the sequence of numbers {c^, a^}™, c;, <Zj = 1 . . . 3. 
The details are given in [Hj. 

Besides the problem related to the color treat- 
ment, the summation over different flavors is also 
a very important problem when the flavor of par- 
tons at the final state is unidentified, as usually. 
In that case a Monte Carlo treatment over flavor 
degrees of freedom has been proposed some time 
ago|21j. showing that the purely gluonic contribu- 
tion falls from 45.7% for 3-jet, to 26.6% for 8-jet 
production j2Tj. 



2.2. On-shell recursive equations 

During the last year much progress has been 
made in the understanding of analytical calcula- 
tions of color amplitudes in perturbative Yang- 
Mills theories. Led by an observation of Witten 
121], Britto, Cachazo and Feng (BCF) have pro- 
posed a new recursion relation for tree amplitudes 
of gluons 3; that naturally arrives at the simplest 
known expressions for those amplitudes in terms 
of Weyl - Van der Waerden spinors, with Maximal 
Helicity Violating vertices as building blocks. Ex- 
plicit calculations have been performed using this 
technique (23JIS]; extensions to amplitudes in- 
volving particles from the electroweak sector [2*5] 
have been pursued and a new approach to one 
loop amplitudes has been proposed |26j employ- 
ing MHV vertices and unitarity arguments as well 
as the use of recursive equations |27I28| . 

The BCF recursion relation features some re- 
markable characteristics, among which the on- 
shell analytic continuation of selected off-shell 
propagators, the analytic continuation of two se- 
lected external momenta in the complex plane 
and a decomposition of a color helicity ampli- 
tude into smaller helicity amplitudes with com- 
plex external momenta that doesn't appear to be 
in direct connection with the decomposition in 
Feynman diagrams. 

For n— gluon amplitude the BCF equation, in 
a diagrammatic representation reads as, 



A{l...n) = 




or in mathematical terms, 

A n (l,2,...,(n-l)-,n + ) - 

n — 3 ^ 

E (^WM, 2,... pr 
i=l h=+,- n ' i 

A n ^i(+p-^,i + l,...,n-2,n-lfj 
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where 

Pn,i =Pn+Pl + ■■■+Pi, 

p 1 1 

Pn,i — Pn.i ~t~ (n— 1 P n 1 \ n] 1 ^ n ' 

Pn-1 = Pn-1 - ( rl _l|P„ ^n-l^n, 
Pn = Pn + ( n -\\p n i \ n ] <Vs-lA n . 

where Ai and Ai are spinors and anti-spinors core- 
sponding to the momentum and (i|P|j] = 

The kinematical operation on the momenta of 
the 1st and nth particles, is called the 'hat' oper- 
ation. 

It can be proven that BCF equation can be 
obtained from the Berends-Giele (or DS) recur- 
sive equations, by making use of the following 
points 

• a special gauge choice, that allows the can- 
celation of all contributions where diagrams 
with the first and the last leg meeting in a 
three-vertex, as well as with the first and 
the last leg meeting in a four-vertex with 
another external leg. These diagrams are 
obviously not-present in the BCF decom- 
position of the amplitude, 

• a set of relations guaranteed by the kine- 
matical transformation (the hat 'opera- 
tion') applied to the chosen momenta, 
that exactly takes care of the apparent 
over-counting of certain Fcynman diagrams 
within the BCF decomposition, 

• and finally a gauge identity, that relates a 
'hatted' contribution arising from a three 
vertex, with the un-hatted three- and four- 
vertex contributions. 

-S— ?_ + -I- 

where 




with V uvpa the QCD four-vertex. 



Although BCF equations are very powerful in 
order to obtain 'analytical' results, their numer- 
ical implementation does not show up any real 
gain as compared to the Berends-Giele (or DS) 
ones. In fact for a moderate number of external 
particles 8 < n < 12 their complexity and there- 
fore their CPU-time consumption is substantially 
larger than that of the Berends-Giele equations 2 . 
Nevertheless their usefulness in computing tree- 
as well as one-loop amplitudes is still an unex- 
plored territory. 

3. HELAC/PHEGAS: results 

HELAC/PHEGAS is a computer package that in- 
corporates the DS approach to compute scatter- 
ing cross section for arbitrary process. In the 
computation of the matrix elements, it includes 
all Standard Model particles and interactions, 
both in Feynman and Unitary gauges. There are 
two options to deal with the colored particles, 
namely the color-connection approach, in which 
all color ordered amplitudes are computed, and 
the color-configuration approach in which the full 
amplitude is given, followed by a Monte-Carlo 
treatment of the color summation. In the phase 
space generation and integration, sector, PHEGAS 
is using a multi-channel approach, each Feyn- 
man graph, identified with a potential generation 
channel, followed by an optimization of the a pri- 
ori weights, entering the calculation of the global 
phase-space density. Moreover for multi-particle 
processes where the number of Feynman graphs 
makes the use of a multi-channel approach impos- 
sible, other phase-space generation methods and 
packages, like HAAG |2U and DURHAM are used. 

HELAC/PHEGAS has been used extensively to 
produce physically relevant results 3 . In the se- 
quel we are going to restrict ourselves to two spe- 
cific examples, in order to reveal its potential for 
physics studies. 

The first example refers to the process p p — > 
t i b b b b. From the physics point of view it con- 
sists the irreducible background of ttHH produc- 
tion, which seems interesting in a high-luminosity 
LHC (SLHC) version, for studying HHH cou- 

2 See also 

3 See for instance 1321331 . 
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pling 34 . From the computational point of view, 
it is a challenging process, and a nice example to 
demonstrate the ability of PHEGAS/HELAC to deal 
with QCD processes in a realistic setup. 

The number of Feynman graphs contributing 
to this process is 1454 (for a gg initial state), 
with 5! color-connection configurations. We have 
used the structure functions and a s from PDFLIB, 
CTEQ-4L (LO). Kinematical decays of t -> bW+ 
have been implemented and the following set of 
cuts has been used: > 20GeV, \r]b\ < 2.5, 
AB, > 0.5. The result for the total cross section 
is 1.053 ± 0.073 (fb) @ LHC energy. 

The second example is the computation of 
Fermion-Loop (FL) contributions in six-fermion 
production processes in e + e~ collisions 35 3 6137) . 
It is an explicit example of the use of DS equations 
to compute one-loop amplitudes. This is achieved 
in a rather straightforward way, by adding to the 
tree-order SM vertices, contributions arising from 
1PI graphs at one loop |38I39| . In such a way us- 
ing FORM we were able to calculate all V1V2V3 and 
V1V2V3V4 vertices at one loop, for arbitrary kine- 
matical configuration (all particles off-shell) and 
for all gauge bosons V — ^,Z 1 W ± . Then, as 
is dictated by the well known quantum-field the- 
oretic argument, the amplitude is just the tree- 
order one with the tree-order vertices replaced by 
the generalized one-loop vertices. 

We take as an example the process e~e + — > 

v^udr^ t + . The number of Feynman Graphs 
contributing to the process is 208, whereas the 
number of DS vertices (the steps needed to com- 
pute the amplitude) is 140. We have used the 
following set of kinematical cuts: Ei,E q > 5GcV 
and mii,m q q > lOGeV, and the result for E = 
500GeV is given by a /ab = 54,96(26) at tree 
order a\jab = 57,31(28) at the one loop, with 
a if -factor given by if/100 = 4.28(2). The MC 
data for this particular run are as follows: MC 
points generated: 1 Million (961792), MC points 
used after cuts: 404842, real time of running: 6 
1/2 hours on a very basic PC. 

4. Outlook 

Recursive equations have been proven to be the 
framework for an efficient matrix clement com- 



putation for arbitrary scattering processes. They 
are the basic ingredients towards the construction 
of an automatized generator including NLO cor- 
rections. The fusion with parton-shower genera- 
tors and the understanding of the working of this 
fusion, will be one of the main tasks in the near 
future. Precision calculations will offer the solid 
basis needed for discoveries in future high-energy 
colliders. 
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